THE PRODUCTION OF STRATIGRAPHICAL 
RANGE-DIAGRAMS BY AUTOMATIC 

METHODS 

by IAN E. PENN 


Abstract. A package of computer programs has been developed on the Institute of Geological Sciences’ IBM 1130 
computer which, after accepting and sorting data submitted on a simple input format, draws stratigraphical range 
and abundance diagrams to any desired scale. The range-diagrams may be of first occurrence; last occurrence; 
taxonomic order; any specified order including sorting within sub-groups. An ornamented lithological plot is 
accurately aligned and a correctly ordered, punctuated and type-set list of fossil names may be produced from com¬ 
panion dictionaries. The resultant diagram is suitable for direct publication. It is estimated that, the programs having 
been installed, the computer method produces diagrams in one-tenth of the time and at one-quarter of the cost (within 
the Institute of Geological Sciences) of manual production. 


In his review of the familiar faunal range-diagram, Bursch categorized the commonly 
used kinds of diagram and quoted Moore’s (1948, p. 324) point that the presentation 
of data, often painstakingly collected, in a single tabular manner may obscure the 
full significance of the results (Bursch 1950, pp. 479-484). This is probably because 
‘the format of the chart itself dramatizes whichever viewpoint is being used’ (McLean 
1972, p. 1). Though highly desirable, the production of several diagrams from the 
same data has never been fully practised perhaps because of the time-consuming 
nature of the work which is in itself not strictly palaeontological. Automatic data- 
processing methods are an obvious solution, for example McLean (1972, p. 5) 
devised a semi-automatic system using relatively cheap card-punch and card-sorting 
equipment. 

By mid-1972 the Institute of Geological Sciences Computer Unit had written 
a package of computer programs in 1130 Fortran for the display of cartographic and 
lithostratigraphic data on the commonly used IBM 1130 computer configuration. 
It was then decided to write programs to produce the conventional faunal range- 
diagrams and link them with part of the larger package which drew ornamented 
lithological sections. At the same time dictionaries were established containing full 
fossil names and the appropriate type-setting instructions so that a copy of the species- 
list for a particular diagram could be automatically punched on paper tape to be fed 
into a phototypesetter which can produce, if required, a correctly type-set and 
punctuated faunal list. This list, together with the diagram, results in the automatic 
production of a faunal range-diagram suitable for publication. 

The major features of this palaeontological program package are here out¬ 
lined (text-fig. 1) while full program listings may be obtained on request. Full details 
of the more generalized central part of the package are in Farmer and Johnson (in 
press). 


[Palaeontology, Vol. 17, Part 3, 1974, pp. 553-563.] 
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text-fig. 1. Flow chart of the various program packages. L DATA, M DATA, and F DATA represent 
the input of lithological mineralogical and fossil data respectively. The various identifiers within rectangles 
represent the various program decks. CUF 11 and CUZ 88 produce only line-printer output; CUZ 84 and 

85 may produce paper tape. 

DATA INPUT 

Data input (text-fig. 2) is short and simple. Lithological data (L DATA in text-fig. 1) 
consists of: the thickness of each successive lithological unit followed by the appro¬ 
priate two-letter code corresponding to the desired lithological ornament punched 
on successive cards. For example, text-fig. 2a shows the entries on successive cards 
of a lithological sequence of (in descending order) 25T2 m of argillaceous limestone; 
4-23 m of mudstone; 0T2 m of oolitic limestone; 1-56 m of clay. The whole entry 
occupies no more than the first nine columns of a standard 80-column punch card. 
At present some twenty-five ornaments are available while a zero entry results in 
blank space which may be ornamented manually if desired. 

The fossil data (F DATA in text-fig. 1) or the mineral data (M DATA) are also 
presented in stratigraphical succession. Thus the first card of each stratigraphical 
sample (text-fig. 2b and c) states the number of species or minerals in the sample 
together with the upper and lower depth limits of the sample; and each subsequent 
card corresponds to each of the species or minerals in these samples. The species or 
minerals are not recorded by name but are given a code number (occupying the first 
four columns of the punch card) which is followed, in the case of species, by a simple 
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(1-4) abundance code corresponding to 'present’, 'fairly common’, 'common’, and 
'very common’ (text-fig. 2b). A zero is entered where the occurrence is doubtful. For 
example, text-fig. 2b shows the entries on successive cards of the first three species 
(code numbers, 4, 3, and 56) of twelve species determined from a depth range of 
17*52-17*92 m, which are estimated as being present, very common and possibly 
occurring respectively. 

The abundance code is the only difference between the fossil and mineral data input 
formats. This is because the mineral 'abundance’ may be in the form of percentages 
(text-fig. 2c) and so one extra column is required on the data card. For example, 
text-fig. 2c shows the entries on successive cards of the three minerals (code numbers 
1, 4, and 5) of a sample at 11-45 m depth which have been found by analysis to con¬ 
stitute 85%, 2%, and 10% of the rock. Percentage fossil data can, of course, be input 
on the M DATA format, while the replacement of abundance codes by equivalent, 
arbitrary percentage values would enable semiquantitative lithological and/or faunal 
data to be displayed alongside quantitative data. It is therefore possible to combine 
lithostratigraphical and biostratigraphical data on the same diagram. It will be seen 
then that the bulk of the data, whether in L DATA or M DATA format, occupies 
no more than the first seven columns of a standard data card. 

Dictionaries giving the full fossil name and corresponding number code are 
permanently stored by accession order and by taxonomic order on punch cards 
(text-fig. 2d) as well as on a magnetic disc. It is simple to insert a new card when a name 
has to be changed or added to the list or to shuffle the deck of cards when the taxo¬ 
nomic order has to be changed. 

A third kind of dictionary entry will be seen between the code number and the 
fossil name (text-fig. 2d). This is a letter code corresponding to convenient groupings, 
e.g. B = brachiopod, BR = rhynchonellacean brachiopod. This code enables the 
computer rapidly to break down large lists of species into meaningful sub-groups 
if desired. 

In addition to the code number and full fossil name, type-setting symbols are 
inserted (text-fig. 2 d) so that the output will be in italic or roman type and have upper 


A lithological data input 
12345678901234567890 
25.12AL 
4.23MU 
0.1 20 L 
1 .56CL 


B fossil data input 

12345678901234567890 
12 17.52 17.92 

4 1 

3 4 
56 0 


C mineralogical data input 
12345678901234567890 

3 11.45 11.45 

1 85 

4 2 

5 1 0 


D dictionary input 


1234567890123456789012345678901234567890123456789012345678901234567890 


224 BR£/@R%HACTORHY NCHIA OBSO LET A ( D% AVIDSON /NON #@S%OWERBY) 

text-fig. 2a-d. Each row of data represents the left-hand side of one 80-column data card and shows 
the entries in their correct columns (as indicated by the italicized numbers) to produce the four different 

kinds of input. 
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MIDDLE JURASSIC MACROFOSSILS IN ACCESSION ORDER MIDDLE JURASSIC MACROFOSSILS IN MEMOIR ORDER 


0 br Kallirhynchia sp. 

1 brKallirhynchia superba S.S. Buckman 

2 br Rhynchonelloidella sp. 

3 brRhynchonelloidella smithi (Davidson) 

4 brRhynchonelloidella smithi crassa Muir-Wood 

5 brRhynchonelloidella wattonensis Muir-Wood 
6btsA vonothyris sp. nov. A 

7bts Tubithyris sp. 

8btsWattonithyris sp. 

9bXsWattonithyris cf. fullonica Muir-Wood 
10bts WattonithyriscX. midfordensis Muir-Wood 
11btlO/77rt/7e//a sp. 

12btlOrnithella bathonica (Rollier) 

13btlOrr?/tf7e//a bathonica bathiensis (Rollier) 
X4bX\Ornithella pupa Muir-Wood 
15btlflug/?e/a sp. 

16btlRugitela cadomensis (Deslongchamps) 

17 aXHolzbergia sp. 

18 atHolzbergia schwandorfensis (Arkell) 

19 aXMorrisiceras sp. 

20 atMorrisiceras comma (S. S. Buckman) 

21 aXMorrisiceras krumbecki Arkell 

22 aXMorrisiceras morrisi (Oppel) 

23 aXMorrisiceras sphaera S. S. Buckman 

24 aXMorrisiceras fornicatum (S.S. Buckman) 

25 aXMorrisiceras skipnum S.S. Buckman 

26 aXMorrisiceras korustes S. S. Buckman 

27 azProcerites sp. 

28 azProceritesprogracilis Cox & Arkell 

29 anProcymatoceras baberi (Morris & Lycett) 

30 aXTulites sp. 

31 atTulites subcontractus (Morris & Lycett) 


PLANTAE 

428 hSolenopora jurassica Nicholson 
427 hSolenopora sp. 

563 h'So/enopora-’sp. 

426 ha\gae [frags] 

430lgnwood [frags] 

137lgnlignite [frags] 

PORIFERA 

423 s'Peronidella pistilliformis’ Lamouroux 

424 s'Cliona'sp. 

425 ssponge [frags] 

ANTHOZOA 
188 scCalamophyllia sp. 

176scChomatoseris hemisphericus (Milne Edwards & Haime) 

177 scChomatoseris orbulites (Lamouroux) 

178 scChomatoseris porpites (Wm. Smith) 

175 scChomatoseris sp. 

181 scCladophyllia babeana (d’Orbigny) 

180 scCladophyllia sp. 

187 scConvexastrea sp. 

195 scDimorpharaea defranciana (Michelin) 

199 scEwardsomeandra vermicularis Milne Edwards & Haime 
198 scEwardsomeandra sp. 

174 scMontlivaltia sp. 

196 scThamnasteria neptuni (d’Orbigny) 

197 scThamnasteria terquemi Milne Edwards & Haime 
194 scThamnasteria sp. 


text-fig. 3. Output of fossil dictionary in accession order (to the left) and in taxonomic order (to the right) 
at 30.9.72. The list has been reduced in scale to the smallest possible size available. 


or lower-case letters. The typesetter is fed by tape produced from the dictionaries 
held on the magnetic disc and produces correctly type-set and punctuated lists (text- 
fig. 3) which it can reduce (or enlarge if desired) so that a large list can be copied at 
a convenient size. The input characters for typesetter instruction have been carefully 
chosen so that in other forms of dictionary output (e.g. the fast lineprinter output 
from the computer or the tape-typewriter print of the paper-tape contents) they will 
not be printed. These other forms of dictionary lists are then much easier to read than 
they would be were the redundant type-setting instructions retained. For example, 
text-fig. 2d shows the dictionary entries for a species (code number 224) of a rhyn- 
chonellacean brachiopod (BR). The typesetting instructions for upper and lower¬ 
case @ and % precede the appropriate letters and the symbols # and / precede the text 
to be printed in roman or italic print. The symbol £ instructs the typesetter to delete 
all the entries to its left if required. 


THE PROGRAMS AND THEIR OUTPUT 

The programs available within the package are labelled (text-fig. 1) by their program 
identifiers. CUC 90 plots a series of lithological sections side by side in specified 
horizontal and vertical relative positions (text-fig. 4) thus providing a skeleton correla¬ 
tion diagram. CUC 70 plots individual lithological sections in the correct vertical 
position relative to plots of stratigraphical range data. Both programs provide 
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text-fig. 4. Skeleton correlation diagram as output by CUC 90 in which the sections are correctly positioned 
stratigraphically and geographically. Labels and tie-lines may be added manually. 


a separate list of those horizons which are too thin to ornament at the chosen scale. 
CUZ 10 and CUF 10 read M DATA and F DATA respectively and store them on 
a magnetic disc for accession by the remaining programs. Thus CUF 11 prints out 
a check-list of the M DATA and F DATA stored by CUZ 10 and CUF 10, and 
CUZ 88 prints out, via the dictionaries, the full taxonomic name of all the different 
minerals or species stored by CUZ 10 and CUF 10. 

CUZ 12 and CUF 12 deduce the stratigraphic range of each mineral or fossil stored 
by CUZ 10 and CUF 10 and plot range and abundance by first occurrence or last 
occurrence (text-fig. 5) or by either of those within designated sub-groups (text-fig. 6). 
It will be seen that the horizontal mark (which is proportional in length to species 






















































































































































































possibly occurring 


558 


PALAEONTOLOGY, VOLUME 17 


CJ\ 

O 

f » i 


O 

i 


SPECIMEN 

CO 

o 

J I I I I 


DIAGRAM 

r\o 

o 

J i I-1-L 


O 

—I— 


o 


J 


I 


o o 
ooo 

C 0 
0 0 0 
0 o 


♦ 

♦ 

I 

I 


Gervillella sp. 

Limatulac f helvetica (Oppel) 
Cercomya undu/ata (J de C. Sowerby' 
Camponectes sp. 

Catmula knorn (Voltz) 

* Grammatodon sp. 

* Protocardiasp. 

. r 

♦ 


'‘IT' 


TT 


3 2. 
<o x 
5 o' 


Praeexogyra acuminata (J. Sowerby) 
Ornithella bathonica (Rollier) 

Or nit he II a sp. 

Inoperna plicata (J. Sowerby) 
serpulids [indet.] 

Entolium corneolum (Young & Bird) 
Eomiodon sp. 

Hybodus sp. 

Bositra buchi (Roemer) 

Cucullaea nuarsensis (Cossmann) 
Oxytomacostatum (Townsend) 
Catmula matisconensis (Lissajous) 


Entolium sp. 

Bositra sp. 

Rhynchonelloidella smith/ (Davidson) 
‘Corbula’sp. |juv.] 

Rhynchonelloidella smithi crassa Muir-Wood 
Rhynchonelloidella wattonensis Muir-Wood 
Rastellum (Arctostree) gregareum (J. Sowerby) 
Rhynchonelloidella sp. 

Liostrea sp. 
crinoids [indet ] 

Camptonectes ngidus (J. Sowerby) 
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Sower by a sp. 

Anisocardia bathensis Cox 
Cucullaea sp. 

Tancredia(lsotancredia)sp [juv ] 
Placunopsis social is Morris & Lycett 
Kallirhynchiasp. 


text-fig. 5. Specimen range-diagram as output by CUZ 85 of last occurrence of Middle Jurassic macro¬ 
fossils from a borehole near Bath. 
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SPECIMEN DIAGRAM 


DEPTH IN METRES 


CARBON A TES 


CLAY MINERALS 


OTHER MINERALS 



text-fig. 6. Specimen range-diagram as output by CUZ 12 showing the grouping of minerals. Calcite 
(shown manually ornamented above the 50% level) was determined quantitatively by wet analysis and this 
was used as an internal standard for calibrating the peak heights of the X-ray diffractograms. The analysis 

was not fully comprehensive. 


abundance) enables every species occurrence to be localized. Where occurrence is 
from a range of depth (e.g. when fossils are recorded from ‘Bed X’) the horizontal 
mark becomes ‘enlarged’ to a box whose height corresponds to the range in depth and 
whose width corresponds to the abundance (e.g. in MK 16 on text-fig. 7). The code 
number of each fossil or mineral is plotted for identification but use of program 
CUZ 85 results in the same plots being produced together with a paper-tape list, in 
the correct order, of the full fossil names obtained from the appropriate dictionary 
corresponding to the code numbers on the plot. This paper-tape when fed through 
the phototypesetter produces a correctly typeset list which may then be affixed to the 
diagram in preparation for publication. CUZ 84 produces a paper-tape in the same 
way as CUZ 85 but plots the range and abundance data in fixed-order formats, e.g. 
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text-fig, 7, Upper Jurassic foraminifera from two boreholes in Central England. The data of borehole 
MK 24, penetrating younger strata, have been arranged by last occurrence as output by CUZ 85. Data 
of borehole MK 16, penetrating older strata, have been arranged as output by CUZ 84 in the order dis¬ 
covered in MK 24. Ammonite subzones have been plotted instead of lithology. Species abundance key as 

in previous diagrams. 
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TEX-r-Fic. 7. Upper Jurassic foraminifera from two boreholes in Central England. The data of borehole 
- u u penetrating younger strata, have been arranged by last occurrence as output by CUZ 85. Data 
of borehole MK 16, penetrating older strata, have been arranged as output by CUZ 84 in the order dis¬ 
covered in MR 24 Ammonite subzones have been plotted instead of lithology. Species abundance key as 

in previous diagrams. 
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